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ABSTRACT 

This  progress  report  outlines  our  research  efforts  on  modeling,  analyses  and  optimization  of  a 
novel  class  of  multi-layered  printed  antennas  for  high  gain  applications.  The  proposed  antennas, 
which  are  based  on  a  multi-layer  printed-circuit  version  of  the  conventional  Yagi  array,  are  very 
attractive  for  applications  that  require  high  gain  antennas  in  a  compact  low-profile  package. 
Presence  of  the  dielectric  layers  not  only  hinders  the  need  for  structural  support  of  the  antenna  but 
also  provides  a  few  more  degrees  of  freedom  for  gain  optimization.  In  addition,  these  antennas  can 
be  made  conformal  to  various  shapes  and  surfaces.  Our  research  during  this  interim  period  has 
resulted  in  three  contributions:  1)  We  have  developed  a  numerical  code  for  efficient  electromagnetic 
modeling  of  these  Yagi-like  structures.  A  novel  feature  of  this  code  is  a  new  semi-analytical 
technique  that  speeds  up  the  evaluation  of  the  corresponding  Green ’s  functions  by  a  factor  of  10  or 
higher.  2)  The  feasibility  of  obtaining  high  gain  from  the  proposed  Yagi-like  arrays  is  investigated 
by  performing  a  detailed  parametric  study  for  structures  with  up  to  5  dielectric  layers.  3)  We  have 
developed  an  electromagnetic  optimization  engine  based  on  Evolutionary  Programming.  This  code 
is  applied  to  optimal  design  of  the  printed  Yagi-like  arrays.  Optimization  is  performed  with  respect 
to  lengths  of  driver  and  director  elements  as  well  as  dielectric  constants  and  thickness  of  layers 
under  different  constraints’  criteria.  It  is  shown  that  a  globally  optimized  three  layer  structures  can 
achieve  a  gain  of  13  dBi  or  higher  without  the  need  for  high  permittivity  dielectric  layers. 
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I.  INTRODUCTION 


This  project  focuses  on  three  research  topics  which  are  of  interest  in  various  antenna 
applications  in  communications  and  radar.  The  topics  under  investigations  are  :  1)  Investigation  of 
a  class  of  “Yagi-like”  multilayer  printed  arrays  for  high  gain  applications,  2)  New  design  concepts 
for  low-profile  dual-band  CP  microstrip  antenna  elements,  and  3)  Performance  Analysis  of  multi¬ 
layer  printed  antenna  structures  on  a  curved  surface. 

A.  Objectives  of  the  Present  Research 

The  objective  in  Topic  1  is  an  in-depth  theoretical  as  well  as  experimental  investigation  of  a 
novel  concept  based  on  a  multi-layer  printed-circuit  version  of  the  conventional  Yagi  array  in  air. 
This  structure  is  very  attractive  for  applications  that  require  high  gain  antennas  in  a  compact  low- 
profile  package.  Presence  of  the  dielectric  layers  not  only  hinders  the  need  for  structural  support  of 
the  antenna  but  also  provides  a  few  more  degrees  of  freedom  for  gain  optimization.  In  Topic  2,  the 
objective  is  to  investigate  a  new  class  of  dual-band  CP  microstrip  antennas,  which  are  proximity- 
fed  in  a  multi-layer  dielectric  substrate.  Due  to  their  electromagnetic  feeding  scheme  these  antennas 
are  expected  to  have  low  spurious  radiation  from  feed  structures  and  thus  well-suited  for  array 
applications  in  GPS  and  satellite  communications.  In  Topic  3,  the  objective  is  to  analyze  the 
radiation  performance  of  printed  antennas  in  a  multi-layered  curved  substrate/ground-plane 
environment.  In  particular,  the  effect  of  cylindrical  curvature  on  the  performances  of  the  antenna 
structures  in  Topics  1  and  2  as  well  as  their  RCS  response  will  be  analyzed. 

The  present  report  details  the  progress  made  on  the  research  Topic  1  during  the  time-period, 
October  1997  -  September  1998.  The  technical  contributions  made  during  this  period  include: 

i)  A  computer  code  for  efficient  electromagnetic  modeling  of  printed  Yagi-like  antennas  in  a 
general  multi-layered  medium  has  been  developed.  A  novel  feature  of  this  code  is  a  new  semi- 
analytical  technique  which  speeds  up  the  evaluation  of  the  corresponding  Green’s  functions  by  a 
factor  of  10  or  higher. 

ii)  Feasibility  of  obtaining  high  gain  from  the  proposed  printed  Yagi-like  dipole  arrays  has  been 
investigated  by  performing  a  detailed  parametric  study  for  structures  with  up  to  5  dielectric  layers. 

iii)  An  electromagnetic  optimization  engine  based  on  Evolutionary  Programming  has  been 
developed.  This  code  has  been  applied  to  optimal  design  of  the  printed  Yagi-like  arrays. 
Optimization  has  been  performed  with  respect  to  lengths  of  driver  and  director  elements  as  well  as 
dielectric  constants  and  thickness  of  layers  under  different  constraints’  criteria. 
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B.  Organization  of  the  Report 

This  progress  report  is  organized  as  follows.  Development  of  an  efficient  Moment  Method 
based  computer  code  for  modeling  of  multi-layer  printed  antennas  in  a  general  magneto-dielectric 
substrate  is  described  in  sections  II- A  and  II-B.  The  application  of  this  code  to  study  the  feasibility 
of  obtaining  high  gain  from  the  proposed  multi-layered  Yagi-like  array  is  presented  in  sections  II- 
C  and  II-D.  Development  of  an  electromagnetic  optimization  engine  based  on  Evolutionary 
Programming  and  its  application  to  optimized  design  of  the  Yagi-like  structure  is  discussed  in 
section  III.  Finally,  future  research  work  planned  for  this  project  is  presented  in  section  IV. 


II.  ANALYSIS  OF  “YAGI-LIKE”  PRINTED  ARRAYS  FOR  HIGH  GAIN 

APPLICATIONS 


A.  Proposed  Structures 

The  proposed  multi-layered  printed  antenna  structures,  shown  in  Figures  1  and  2,  present  a 
low-profile  alternative  to  the  conventional  Yagi-Uda  array  of  wire  dipoles  which  are  widely  used  in 
communications  and  radar  for  high-gain  applications.  Although  the  Yagi-array  has  a  relatively 
simple  structure  made  of  thin-wire  elements,  its  overall  length  could  become  very  large  for  high 
gain  and/or  low-frequency  applications  mainly  due  to  the  large  number  of  directors  and  the 
required  spacing  between  the  elements.  In  addition,  a  mast  must  be  used  to  structurally  support  the 
array  elements  in  air.  The  structure  in  Figures  1  is  a  multi-layer  printed-circuit  version  of  the  Yagi 
array  of  wire  dipoles  in  air.  This  structure  consists  of  a  reflector,  a  driver  and  a  finite  number  of 
embedded  director  strip  elements.  An  alternative  structure,  shown  in  Figure  2,  is  a  microstrip 
version  of  the  array  in  Figure  1  and  has  a  ground-plane  in  place  of  the  reflector  element.  These 
structures  give  more  degrees  of  freedom  (i.e.,  the  dielectric  constants  of  the  layers)  when 
optimizing  their  design.  They  can  also  be  made  conformal  to  various  shapes  and  surfaces  (e.g., 
missiles,  aircraft,  or  other  flying  objects).  In  addition,  unlike  the  gain-enhancement  techniques 
reported  for  an  embedded  microstrip  dipole  [1,2],  the  proposed  structures  can  achieve  high  gains 
without  the  need  for  very  high  permittivity  dielectric  layers. 

The  above  concepts  can  also  be  extended  to  other  patch  geometry  (rectangular,  circular,  etc.) 
and  feed-mechanism  (coaxial,  microstrip-line,  etc.)  as  shown  in  Figure  3.  It  may  also  be  possible 
to  design  these  “Yagi-Like”  structures  for  applications  (e.g.,  satcom)  that  require  circular- 
polarization  (CP),  by  using  for  the  driver  and  director  elements  such  structures  as  printed  cross¬ 
dipoles,  cross-fed  patches,  nearly  square  patches,  slotted  patches,  etc.  Finally,  a  phased  array 
consisting  of  the  proposed  Yagi-like  antenna  elements  may  be  formed  for  beam-steering  purposes 
as  depicted  in  Figure  4. 
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Figure  1:  A  printed  Yagi-like  dipole  Array 
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Figure  2:  A  Yagi  stacked  microstrip  array  in  a 
multi-layered  substrate 


Figure  3:  A  coaxial-fed  Yagi-like  stacked 
microstrip  patch  array 


Figure  4:  A  phased-array  of  Yagi-like  stacked 
patch  sub-array  elements 


B.  Development  of  an  Electromagnetic  Model  for  the  Antenna  Structures 

To  optimize  and  design  the  proposed  Yagi  structures,  an  efficient  numerical  engine  for 
accurate  modeling  of  a  general  multi-layer  printed  antenna  structure  is  needed.  For  the  present 
problem  the  CPU  efficiency  of  a  numerical  engine  is  of  particular  importance  since  one  has  to 
optimize  the  antenna  gain  with  respect  to  the  reflector,  driver  and  directors  lengths,  thickness  of  the 
layers  and  their  dielectric  constants.  A  typical  multi-parameter  global  optimization  technique,  may 
require  the  engine  to  be  run  hundreds  of  times  in  order  to  achieve  a  desired  gain  specification. 

In  this  work,  we  have  developed  an  efficient  moment-method  based  technique  that  is  well- 
suited  for  the  present  problem.  The  technique  uses  the  so-called  mixed-potential  integral  equation 
(MPEE)[3,4]  which  in  general  can  model  planar  microstrip  structure  of  arbitrary  shape  in  a  multi¬ 
layered  medium  (Figure  5).  Our  implementation  of  MPDE  utilizes  a  newly  developed  technique  for 
semi-analytical  evaluation  of  the  corresponding  Green’s  functions  and  is,  therefore,  well  suited  for 
a  medium  with  many  dielectric  layers. 
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Figure  5:  A  multi-layered  microstrip  antenna  structure 


In  general,  the  MPIE  for  a  multi-layered  structure  can  be  written  as: 

N 


i  =  1  si 


(*■*)  JJ*)  +  TiVtG(EP! (*■*■)  VfJJ*  ) 

K-0 


dS’i 


-  ZspJsp(x)  =  -EZ> 


(1) 


x  -  (x,y)  e  S  ;  p  -  1,2,....  N 


where  Jsj  and  Zsi  are  the  electric  surface  current  and  the  surface  impedance  of  Sj,  respectively; 


Einc^P)  is  the  tangential  (incident)  electric  field  impressed  at  the  p-th  surface.  G^.  and  are  the 
Green's  Functions  of  Magnetic  and  Electric  types,  respectively,  evaluated  at  the  p-th  layer  due  to 
a  source  at  the  i-th  layer.  The  Green's  functions  are  of  the  Sommerfeld  integral  type  which  can  be 
expressed  as  the  inverse  Bessel  transformation, 


GJz,P)  =  2k  \G  Jz,X)J0(k0pX)XdX  ;  p  =  ^/(x-x')2  +  (y-y')2  (2) 

0 

(1) 

where  the  spectral-domain  Green's  functions  Gu(z,X)  satisfy  the  voltage  distribution  in  an 
equivalent  multi-section  transmission  line. 


B.l.  A  Semi-Analytical  Method  for  the  Evaluation  of  Green’s  Functions 

In  a  typical  moment  method  solution  of  the  MPIE  in  (1),  one  has  to  numerically  compute  the 
spatial-domain  Green's  functions  GE  and  GM  in  (2).  The  integral  in  (2),  however,  is  of  the 
Sommerfeld  type  which  is  not  trivial  to  integrate  due  to  its  surface  wave  and  branch-cut 
singularities.  For  an  N  layered  structure  with  metalizations  on  all  layers,  there  are 

N 

M  =  2l{N-m  +  i)  unique  Green's  functions,  requiring  M  numerical  evaluations  of  the  Sommerfeld 

m-1 
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integral  for  each  kQp  value,  i.e.,  the  normalized  lateral  distance  between  a  source  and  an 

observation  point.  Each  spectral-domain  G  has  to  be  recursively  determined  from  an  equivalent 
multi-sectional  transmission  line  representing  the  multi-layered  medium.  Since  a  large  number  of 
kQp  values  are  required,  the  CPU  time  spent  on  computation  of  the  Green's  functions  can  be 

prohibitively  large  when  one  attempts  to  develop  a  simulation  /  optimization  program  for  the  design 
of  microstrip  arrays  located  in  such  an  environment. 

In  this  work  we  have  developed  a  semi-analytical  method  for  evaluation  of  the  Green’s 
functions  which  significantly  speeds  up  the  numerical  evaluation  of  these  integrals.  Here  we  briefly 
present  an  outline  of  the  method.  A  more  detailed  formulation  is  presented  in  Appendix  A.  In  this 
semi-analytical  method  we  first  deform  a  portion  the  contour  of  integration  above  the  real-axis  into 
the  complex  X  .  plane,  in  order  to  avoid  the  singular  behavior  of  the  integrand  near  any  surface 
wave  poles,  and  then  back  to  the  real  axis  to  infinity.  Cubic-spline  interpolations  are  used  to  curve- 
fit  the  spectral-domain  integrand,  G ,  in  each  region,  except  in  the  tail  of  the  integrand,  i.e.,  Xmax 
(  »1  )  to  infinity  region,  for  which  the  integration  range  is  divided  into  small  pieces,  using  an 
adaptive  marching  scheme,  and  G  is  approximated  on  each  piece  by  the  first  few  terms  in  it's 

Taylor-series  expansion  in  X *2.  The  resulting  integrals  can  then  be  expressed  in  closed  forms  in 
terms  of  known  special  functions,  namely,  Bessel  and  Struve  functions.  These  functions  are  of 
distance  p  only  and  independent  of  the  dielectric  layers'  thickness  and  material  properties. 
Therefore,  once  the  constants  in  the  cubic-spline  interpolations  and  Taylor-series  expansions  are 
computed  and  stored,  both  the  electric  and  magnetic  type  Green’s  functions  can  be  evaluated  for  all 
values  of  p  at  once,  resulting  in  a  substantial  reduction  in  the  overall  computation  time.  In 
addition,  a  scheme  is  developed  in  which  the  resulting  special  functions  are  computed  only  once 
for  a  given  set  of  kop  values,  stored  and  then  used  for  all  M  Green’s  functions  of  an  N-layer 

structure.  This  results  in  an  additional  CPU  time  saving  when  metalized  layers  are  present  on  two 
or  more  of  the  dielectric  layers. 

As  an  example,  the  Green’s  function  Ge  and  Gm  of  a  five  layer  structure  are  compared  with 
the  exact  numerical  integration  in  Figure  6.  The  results  are  shown  for  the  case  when  the  source  and 
observation  points  are  located  on  the  2nd  and  3rd  layers,  respectively,  The  integral  averaging 
technique  [5]  is  used  to  speed  up  the  numerical  integration  for  the  exact  results.  Excellent 
agreement  for  both  the  real-axis  as  well  as  the  deformed  complex  parts  of  the  contour  are  obtained. 
The  normalized  CPU  time  required  for  calculating  all  30  distinct  Green’s  functions  of  this  five 
layer  structure  at  100  kop  points  was  only  10  seconds  using  the  present  method  as  compared  to 

more  than  220  seconds  when  the  numerical  integration  was  used.  In  general,  as  shown  in  Figure 
7,  the  CPU  time  saving  increases  as  the  numbers  of  dielectric  layers  and  the  observation  point 
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Magnitude 


increase,  indicating  the  efficiency  of  this  method  for  antenna  structures  with  a  large  number  of 
conductors  and  dielectric  layers. 


Kop 

Figure  6:  Results  for  the  real-axis  portion  of 
the  contour.  erI  =  2,  er2=  3,  er3=  4,  er4=  6,  er5  = 
2j  hodj  =  0.1,  k0d2  =  0.1,  kod3  0.05,  kod4 
0.1,  Ms  =  0.2.. 


Number  of  KoP  Points 

Figure  7:  CPU  time  comparison  between 
numerical  and  semi-analytical  methods 


B. 2.  Moment  Method  Implementation 

To  solve  the  integral  equation  (1)  numerically  for  the  unknown  current  densities,  Jsi, 
i=l,2,...N,  we  have  used  a  Galerkin/Moment  Method  solution.  For  the  configurations  in  Figures  1 
and  2  each  strip  conductor  is  assumed  to  have  an  electrically  small  width,  thus  the  current  on  the 
strip  is  primarily  flowing  along  its  length,  i.e.,  in  the  x-direction.  Therefore,  each  dipole  is  divided 
into  20  rectangular  cells  per  wavelength  in  x-direction  and  only  1  cell  in  y-direction.  On  each  cell  a 
constant  basis  function  is  used  for  the  current.  Constant  basis  functions  are  also  used  for  the 
charge  distribution,  VfJs-  The  charge  cells,  however,  are  shifted  by  1/2  cell  with  respect  to  the 
current  cells  in  order  to  satisfy  the  continuity  equation  in  an  average  sense. 

C.  High  Gain  Conditions  for  a  Printed  Dipole  in  a  Multi-Layer  Medium 

In  order  to  understand  the  physics  of  radiation  from  printed  Yagi-like  antennas  in  Figures  1 
and  2,  we  first  need  to  examine  the  conditions  of  obtaining  high  gain  from  a  printed  dipole  in  the 
absence  of  any  director  elements.  Let  us  consider  a  Hertzian  dipole  embedded  in  a  multi-layered 
magneto-dielectric  medium  (see  Figure  8).  There  are  two  types  of  high  gain  conditions  for  this 
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structure  [2].  Type  I  condition  requires  that  odd  layers  have  high  wave-impedance  (i.e.,  large  /lr ) 
,  even  layers  have  low-wave  impedance  (i.e.,  large  er)  and 


Tljd-I  _  1  niZo  _  1  .  nAi  _  ^  -•  o 

-  -  )  «  >  *  >  *  -  ** 

\0  2  \0  4  A0  4 


(3) 


Type  II  condition  requires  that  odd  layers  have  low  wave-impedance  (i.e.,  large  er)  ,  even  layers 
have  high  wave-impedance  (i.e.,  large  fir)  and 


nIdl  _  1  rtjZp  _  1 


^-  =  4  .  i>2 


(4) 


Figure  8:  An  embedded  dipole  in  a  multi-layered  medium 

where  n,  =  is  the  refractive  index  of  layer  i.  It  is  obvious  that  for  a  purely  dielectric 

medium  (i.e.  /J,=l ) ,  the  dielectric  layers  should  alternate  between  low  and  high  permittivities.  For 
example  as  shown  in  [6],  for  a  resonant  microstrip  dipole  in  a  two-layer  medium  with  erl  =  2,  type 
I  condition  requires  an  er2  of  larger  than  10  if  a  high  gain  of  larger  than  10  dBi  together  with  high 
radiation  efficiency  (i.e.,  minimum  of  power  lost  to  the  surface-waves)  is  desired.  This  may  not  be 
feasible  in  practice  because  of  the  cost  and  heavy  weight  of  high  permittivity  dielectric  materials. 
As  will  be  shown  in  the  next  section,  introduction  of  director  elements  on  the  interfaces  between 
dielectric  layers  resolves  this  problem. 

D.  Results  and  Discussions 

A  three  layer  “Yagi-like”  stacked  microstrip  dipole  array  is  shown  in  Figure  9.  Figure  10 
shows  the  gain  of  this  structure  plotted  versus  the  director  length  (  L2  =  L3 )  in  X0,  the  free-space 
wavelength.  All  three  layers  are  assumed  to  have  £rl  =  er2  =  £r3  and  equal  thickness,  d,=  d2=  d3  = 
0.25  Xd,  where  A,d  is  the  wavelength  in  dielectric.  Width  of  each  strip  element  is  W  =  0.021  X0  and 
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the  length  of  the  driver  is  optimized  for  a  resonance  at  the  normalized  frequency  of  f  /  f0=l.  As  can 
be  seen  for  er  =  2  a  gain  of  about  9.6  dBi  is  obtained  when  the  each  director  has  a  length  of  about 
0.32  free-space  wavelength ;  the  length  of  the  driver  in  this  case  is  Lt  =  0.302Xq  .  We  note  that 
the  same  structure  without  the  directors  but  with  identical  dielectric  parameters  has  a  gain  of  only 


about  5.84  dBi  [6]. 


Figure  9:  A  three  layer  microstrip  Yagi-like 
antenna 


Figure  10:  Gain  of  a  printed  3-element 
stacked  microstrip  Yagi;  erl  =  =  £,.3  =2  and 

4  \  d,=  d2=  d3==0.25  Xd 


The  Yagi-like  structure  can  be  further  optimized  and  made  more  compact  by  keeping  d,  fixed 
and  varying  d2=  d3,  as  shown  in  Figure  11.  As  seen,  for  d2=  d3  =  0.1  Xd  and  X2=  X3=  0.36  X0 ,  a 
gain  of  about  10.1  dBi  is  obtained  with  an  overall  thickness  of  0.45  Xd  (nearly  0.32  X0  ).  The 
corresponding  driver’s  length  is  L,=  0.296  X0.  The  radiation  efficiency  (i.e.,  ratio  of  the  radiated 
power  to  the  total  of  the  radiated  plus  the  surface-wave  powers)  in  this  case  is  about  90%  as 
compare  to  only  about  67%  for  an  embedded  dipole  case  without  the  director  elements.  The  Yagi- 
effect  is  in  fact  at  work:  addition  of  the  two  directors  tends  to  focus  the  energy  in  the  broadside 
direction  and  minimize  the  power  lost  to  the  surface-waves.  Another  evidence  of  this  is  the  plot  of 
the  current  distributions  in  Figure  12  where  as  seen  a  relatively  substantial  amount  of  current  is 
induced  on  the  director  elements.  The  trade-off  in  this  case  is  the  gain  bandwidth,  which  as  shown 
in  Figure  13,  is  much  narrower  for  the  Yagi-like  structure  than  for  the  embedded  dipole  case 
without  the  director  elements.  In  Figure  14,  the  gain  versus  the  director  length  is  shown  for  a 
structure  with  d!  =  0.25  Xdl  ,  er]  =  er2  =  2  and  er3  =  4  and  for  various  values  of  d2/  Xd2  =  d3/  Xd3. 
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Again,  a  high  gain  of  more  than  1 1  dBi  is  obtained  when  d2/  Xd2=  d3/  Xd3=  0.235  and  L J  Xd2  =  L3/ 
Xtf  =  0.36.  The  corresponding  resonant  length  of  the  driver  in  this  case  is  about  O.288A,0.  The 
same  structure  without  the  directors  has  a  gain  of  only  about  8  dBi.  For  this  example  the  array  gain 
decreases  as  values  of  d2  and  d3  are  further  increased. 


0.05  0.15  0.25  0.35  0.45 

Director  Length  /  Wavelength 

Figure  11:  Gain  of  a  printed  3 -element  stacked 
microstrip  Yagi;  £ri  =  er2  =  £r3  -2,  d]  =  0.25  A,d 
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Figure  13:  Gain  of  a  3-element  stacked  microstrip 
Yagi  array  as  a  function  of  normalized  frequency, 
dj  =  0.25  Xd 
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Figure  12:  Current  distributions  on  the  stacked 
Yagi  elements;  =  er2  =  er3  =2,  dj  =  0.25  Xd 
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Figure  14:  Gain  of  a  3 -element  stacked 
microstrip  Yagi  as  a  function  of  L2/  Xd2  = 
L3/  X.d3;  The  label  on  each  curve  is  the 
ratio:  d2/  Xd2  =  d3/  Xi3. 


Higher  gain  values  can  be  obtained  by  increasing  the  number  of  dielectric  layers  and  combine 
the  high  gain  condition  discussed  in  section  H.C  (i.e.,  alternate  between  low  and  high  £r  )with  the 
Yagi  effect  (i.e.,  directors  at  each  interface).  Figure  15  shows,  the  results  for  a  five  layer  antenna 
with  d,  =  0.25  A,dl ,  erl  =  £r2  =  £r4  =  2,  £r3  =  er5  =  4  and  for  various  values  of  d2/  Xd2  =  d3/  Xd3= 
d J  X,d4=  d5 /  X.d5.  Again,  we  have  assumed  all  directors  are  of  equal  length.  As  can  be  seen,  a  high 
gain  of  more  than  14  dBi  is  obtained  when  d2/  Xd2  =  d3/  Xd3  =  d4/  Xd4  =  d5/  Xd5  =  0.245.  An 
inspection  of  the  current  distributions  for  this  case  shows  that  the  magnitude  of  the  current  on 
directors  1  and  3  are  substantially  more  than  those  induced  on  the  directors  2  and  4.  It  is 
noteworthy  that  directors  1  and  3  are  at  the  transition  interface  between  a  low  and  a  high 
permittivity  layers,  while  directors  2  and  4  are  at  the  interface  between  a  high  and  a  low  permittivity 
layers.  This  phenomenon  should  be  explored  in  more  details  in  the  future  phases  of  this  research 
since  it  may  provide  some  guidelines  for  optimal  design  of  these  high  gain  antennas. 


Figure  15:  A  five-layer  microstrip  Yagi-like 
antenna.  Gain  is  plotted  as  a  function  of 
director  length;  the  label  on  each  curve  is  the 
ratio:  d2/  =  d3/  Xd3.  =  d4/  X43—  d$/  . 

Current  distributions  are  plotted  for  the 
maximum  gain  case  with  d2/  Xd2  =  d3/  A,d3 
—  d4/  A.43=  X,(]5 =  0.245  and  L2/  Xd2  ~  L 3/  ^d3 

“  L4/  A.d4  =  L5 /  =  0.355. 


Director  Length  /  Wavelength  in  Dielectric 
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In  all  of  the  above  examples  the  thickness  of  the  first  layer  was  assumed  to  be  d ,/  Xdl  =  0.25 
The  feeding  mechanism  in  many  microstrip  patch  antenna  applications  is  coaxial  or  probe-fed 
which  makes  it  undesirable  to  have  a  thick  substrate  thickness  due  to  spurious  radiation  from  the 
probe.  As  we  have  shown  in  Figure  16,  however,  it  is  possible  to  obtain  a  high  gain  even  with  an 
electrically  thin  substrate.  For  this  case  the  frequency  is  1  GHz  and  the  substrate  thickness  is  fixed 
at  dj=  2.5  mm.  A  gain  of  more  than  11  dBi  is  obtained  when  d2/  Xd2  =  d3/  Xd3  =  0.23. 
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Figure  16:  A  three-layer  microstrip  Yagi  with  a  thin  substrate.  Gain  is  plotted  as  a  function  of 
director  length;  the  label  on  each  curve  is  the  ratio:  d2/  Xd2  =  d3/  Xd3  . 
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III.  GAIN  OPTIMIZATION  USING  EVOLUTIONARY  PROGRAMMING 


The  attempts  outlined  in  section  1I-D  to  optimize  the  Yagi-like  structure  for  high  gain  were 
arbitrary  in  nature  and  were  done  using  a  parametric  study.  Nevertheless,  it  was  shown  that  a  three 
layer  structure  with  er  =  2  can  produce  a  gain  as  high  as  10  dBi  or  more,  indicating  its  potential  for 
high  gains.  Further,  our  investigation  has  indicated  that  the  objective  function  relating  the  gain  to 
the  antenna  parameters  is  highly  non-linear  and  possesses  multiple  local  optima.  To  successfully 
optimize  this  structure  multi-parameter  global  optimization  schemes,  such  as  evolutionary 
algorithms,  are  necessary. 

A.  Evolutionary  Programming 

The  evolutionary  algorithms,  such  as  genetic  algorithms  (GAs)[7],  evolutionary 
programming  (EP)[8],  and  evolutionary  strategies  (ES),  are  heuristic  population-based  search 
procedures  that  incorporate  random  variation  and  selection.  Of  the  three  paradigms  GAs  are  well 
known  to  the  electromagnetics  community  (see  [9]  for  a  detailed  bibliography).  Even  though 
several  successful  applications  have  been  reported,  recent  research  has  identified  some 
inefficiencies  in  GA  performance.  This  degradation  in  efficiency  is  apparent  in  applications  with 
highly  epistatic  objective  functions,  i.e.,  where  the  parameters  being  optimized  are  highly 
correlated.  On  the  other  hand,  EP  and  ES  are  more  robust  to  epistatic  objective  functions  and 
coordinate  rotations.  EP  has  been  shown  to  be  more  efficient  than  GA  on  many  function 
optimization  problems  [8]. 

The  dynamics  of  GA  have  been  explained  through  the  building  block  hypothesis  and  the 
schema  theorem  [7],  which  are  not  fully  accepted  in  the  evolutionary  computation  literature.  On  the 
other  hand,  the  convergence  theory  for  EP  is  well  established  and  EP  has  been  proven  to 
asymptotically  converge  to  the  global  optimum  with  probability  one,  under  elitist  selection  [8]. 
Further,  EP  is  well  suited  for  optimizing  continuous,  discrete,  and  mixed  parameter  optimization 
problems.  Binary  GAs  require  the  parameters  to  be  coded  as  bits;  this  could  limit  the  accuracy  of 
the  solutions  obtained.  Further,  the  selection  of  the  crossover  and  mutation  probabilities  is  quite 
arbitrary  and  they  are  not  adapted  during  evolution.  The  selection  of  the  initial  values  for  the 
strategy  parameters  for  EP  and  ES  are  well  defined  and  efficient  adaptive  and  self-adaptive 
techniques  exist  for  adapting  these  parameters  during  evolution. 

The  algorithm  for  a  standard  EP  consists  of  seven  general  steps: 

1 .  Initialization:  Initialize  parents  uniformly  over  the  population  size  of  N. 

2.  Fitness  Evaluation  of  Parents:  evaluate  the  objective  function  for  each  parent. 

3.  Mutation:  From  each  parent  a  single  offspring  is  generated  by  adding  a  Gaussian  random 
variable  with  the  standard  deviations,  SD. 
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4.  Fitness  Evaluation  of  Offspring:  evaluate  the  objective  function  for  each  offspring;  The  parents 
and  offsprings  constitute  the  new  population. 

5.  Tournament:  each  solution  competes  against  M  ( <  N )  randomly  chosen  opponents  and  the 
number  of  wins  for  each  solution  is  obtained. 

6.  Selection:  the  solutions  with  the  greatest  number  of  wins  are  selected  to  be  the  parents  of  the 
next  generation. 

7.  Termination:  is  the  best  parent  satisfactory?  If  yes,  terminate,  otherwise,  jump  to  step  2. 

A  computer  code  in  C++  was  developed  based  on  the  above  procedure.  Before  applying  this  code 
to  the  gain  optimization  of  the  multi-layered  printed  array,  the  accuracy  of  the  code  was  first 
verified  by  applying  it  to  the  optimization  of  three  well-known  and  challenging  antenna  problems: 
thinned  phased  array,  Yagi-Uda  array  and  aperiodic  array  [10,11].  These  are,  respectively, 
discrete,  continuous  and  mixed  parameter  optimization  problems.  To  the  best  knowledge  of  the  PI, 
these  efforts  are  the  first  applications  of  evolutionary  programming  in  electromagnetics. 

In  the  thinned  phased  array  problem,  the  objective  is  to  turn  off  some  elements  in  a  uniformly 
spaced  or  periodic  array  (with  an  inter-element  spacing  of  one-half  wavelength)  to  create  a  desired 
side-lobe  level  (SLL).  For  a  2N-element  symmetric  array,  thinning  has  22N  possible  combinations 
which  is  substantially  large  for  large  N.  We  applied  EP  to  optimize  a  200  element  array  for  a  SLL 
of  -20  dB.  The  solution  obtained  for  the  200  element  array  had  124  elements  turned  on  and  an  SLL 
of -20.21  dB.  The  best  solution  in  the  Genetic  Algorithm  (GA)  literature  had  145  elements  turned 
on  [12],  which  has  17  percent  more  number  of  elements.  For  the  case  of  a  six  element  Yagi-Uda 
array,  EP  optimization  with  respect  to  the  elements’  lengths  and  inter-element  spacings  resulted  in  a 
directivity  of  13.70  dBi  which  is  comparable  to  the  directivity  of  about  13  dBi  reported  in  the 
literature  using  conventional  optimization  techniques. 

B.  Implementation  of  EP  for  the  Printed  Yagi-Like  Structure 

To  optimize  the  gain,  the  stacked  printed  dipole  array  is  represented  as  a  vector  consisting  of 
the  length  of  the  driver  element,  the  director  lengths,  the  dielectric  constants,  and  the  dielectric 
thickness  as  given  below: 

X  =  [xref,  xdri,  xjl) ...  xdir(N-l),  epr(l) ...  epr(N),  dol(l) ...  dol(N)  ]T  (5) 

where  xref  and  xdri  are  the  lengths  of  the  reflector  and  driver,  and  xdir( i)  ,  epr( i)  and  dol( i)  are  the 
length  of  the  i-th  director  element  in  free-space  wavelength,  7,0,  the  dielectric  constant  of  i-th  layer 
and  the  thickness  of  the  i-th  layer  in  dielectric  wavelength,  7^,  respectively.  In  order  to  ensure  that 
the  optimized  solutions  obtained  were  practically  feasible,  the  following  constraints  were  chosen  : 
a)  Input  resistance,  Rin,  should  be  larger  than  a  minimum  acceptable  value,  b)  First-layer  dielectric 
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thickness  should  be  larger  than  a  minimum  value;  for  the  microstrip  version  of  the  structure,  a  very 
thin  substrate  results  in  a  very  small  bandwidth  (alternatively,  a  requirement  on  the  input 
impedance  and/or  gain  bandwidths  can  be  substituted  for  this  constraint),  c)  Total  thickness  of  the 
structure  should  be  smaller  than  a  maximum  value  in  order  to  ensure  a  more  compact  and  low 
profile  design.  The  constraints  were  implemented  by  penalizing  solutions  that  violated  these 
constraints.  The  objective  function  (which  is  to  be  minimized)  is  then  given  by 

Fitness  =  -  GaindB  +  Pa  +  Pb  +  Pc  (6) 

where  GaindB  is  the  gain  in  dBi  obtained  from  the  moment-method  solution  of  the  current 
distributions  in  MPIE  in  (1)  and  Pa ,  Pb  and  Pc  are  the  penalty  criterion  for  violating  the  above 

constraints  in  (a),  (b)  and  (c),  respectively.  In  addition,  limits  were  set  on  the  maximum  allowable 
values  of  each  layer’s  dielectric  constant  and  thickness  during  the  optimization  process  as 
discussed  below. 

In  the  optimization  process  the  driver  length  was  set  to  its  the  resonance  length  (using  a 
conventional  quasi-Newton  optimization  scheme)  in  order  to  make  the  input  impedance,  Zin,  purely 
resistive,  where  as  the  xref ,  xdir(i)  ,  epr(i)  and  dol(i)  were  optimized  using  EP.  A  standard  EP 
algorithm  for  continuous  parameter  optimization  was  applied  as  depicted  in  the  flow-chart  shown 
in  Figure  17. 

C.  Results  and  Discussions 

A  number  of  simulations  for  a  3-layer  version  of  the  microstrip  structure  in  Figure  2  were 
performed.  Four  cases  are  reported  here.  For  the  first  three  cases  dielectric  constants  were  set  fixed 
to  £rl  =  er2=  £  r3=  2  with  constraints  of  minimum  of  0.1  Xd  for  dt  and  0.05  Xd  for  d2and  d3.  Width 
of  each  strip  element  was  w  =  0.02X0  .  For  Case  1  the  solutions  were  optimized  without  any 

constraints  on  Rin  and  the  total  thickness  of  the  structure.  High  gains  of  11  to  13  dBi  and  higher 
were  quickly  obtained  within  35  generations  (350  fitness  evaluations).  Table  1  shows  the 
optimized  parameters  for  a  sample  result  obtained  after  only  10  generations.  For  this  case  a  gain  of 
13.57  dBi  is  obtained  with  a  radiation  efficiency  of  93%,  which  means  a  minimum  amount  of 
power  is  lost  to  the  surf  ace- waves.  Due  to  the  lack  of  constraint  on  the  input  resistance,  however,  a 
very  low  value  of  Rin  is  obtained.  The  current  distributions  for  this  case  are  plotted  in  Figure  18. 
The  director  1  plays  a  major  role  in  enhancing  the  gain  as  is  evident  from  the  magnitude  of  the 
current  induced  on  its  surface.  Nevertheless  as  shown  in  table  2,  presence  of  both  directors  is 
crucial  in  obtaining  the  optimum  gain,  an  indication  that  the  Yagi  effect  is  at  work- 
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Output 


Figure  17:  Implementation  of  the  EP  algorithm  for  optimization  of  the  printed  Yagi-like  array. 


18 


Solutions  in  Case  2,  were  obtained  subject  to  an  additional  constraint  of  Rin  >10  ohms.  This 
constraint  causes  a  drastic  increase  in  the  number  of  generations.  Gains  of  10  to  11.5  dBi  were 
obtained  within  about  140  generations.  Optimized  parameters  resulted  from  the  EP  optimization 
after  about  32  generations  are  listed  in  Table  1.  Finally  in  Case  3,  in  addition  to  the  above 
constraints,  a  constraint  was  also  put  on  the  maximum  total  thickness:  D  =  dj  +  d2  +  d3  <  0.5  Xd. 
Within  about  150  generations  a  high  gain  of  about  9.88  dBi,  with  D  =  0.447  A,d  (nearly  0.316  X0) 

was  obtained.  The  corresponding  current  distributions  are  plotted  in  Figure  19.  The  constraint  on 
Rin,  forces  the  magnitude  of  the  current  on  the  feeder  element  to  be  less  than  those  induced  on  the 
director  elements.  Removal  of  both  directors  in  this  case  results  in  again  of  only  about  4.30  dBi. 


Table  1:  Optimized  3-layer  ‘Yagi-like’  stacked  microstrip  array 


Gain 

(dBi) 

dl 

(K) 

d2 

(^d) 

d3 

(^•d) 

xdri 

(x.) 

xdri(l) 

On) 

xdri(2) 

(K) 

1 

BjjE^3EM19 

1 

Case  1 

13.57 

0.240 

0.120 

0.360 

0.3186 

0.3180 

0.3889 

0.109 

93% 

Case  2 

11.10 

0.299 

0.252 

0.491 

0.2955 

0.3005 

0.3327 

14.43 

96% 

Case  3 

9.88 

0.125 

0.092 

0.230 

0.2838 

0.2881 

0.3686 

20.32 

81  % 

Case  4 

15.66 

0.192 

0.299 

0.258 

0.3198 

0.1583 

0.4692 

30.18 

97% 

Table  2:  Effects  of  director  elements  on  gain  in  Case  1 


Both  Directors 
present 

Director  1 
removed 

Director  2 
removed 

Both  Directors 
removed 

Gain 

(dBi) 

13.57 

2.35 

8.64 

4.30 

Figure  18:  Current  distributions  for  Case  1 


Length  (cell  #) 

Figure  19:  Current  distributions  for  Case  3 
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Behavior  of  these  three  cases  as  a  function  of  normalized  frequency  is  shown  in  Figure  20. 
As  can  been  seen,  Case  1  with  no  constraints  on  Rin  and  D  results  in  the  highest  gain  at  the 
expense  of  a  very  low  gain  bandwidth.  Case  2  and  3  with  a  constraint  on  Rin,  however,  result  in 
much  higher  bandwidths. 


Figure  20:  Gain  of  a  printed  3-elerfient  stacked  microstrip  Yagi;  £rl  =  £r2  =  £r3  =  2. 

Finally  in  Case  4  the  dielectric  constants  were  allowed  to  vary  in  the  range  of  1.5  to  13  with  a 
constraint  of  Rin  >10  ohms.  It  was  noted  that  £rl  and  £r2  always  tended  to  1.5  while  £r3  tended  to 
13;  this  is  consistent  with  the  high  gain  condition  for  an  embedded  dipole.  A  gain  of  15.66  dBi 
and  radiation  efficiency  of  97%  were  obtained  with  £rl  =  £r2  =1.5  and  £r3  =  11.23  after  114 
generations.  The  corresponding  optimized  parameters  are  shown  in  Table  1. 


IV.  FUTURE  WORKS 

Research  works  on  the  printed  Yagi-like  antennas  that  are  either  presently  in  progress  or  are 
planned  for  the  1998-1999  academic  year  include: 

-  Extension  of  the  moment-method  analysis  of  section  II-B  to  the  Yagi-like  microstrip  patch 
array,  fed  with  a  coaxial  line  ( see  Figure  3). 

-  Theoretical  investigation  of  Gain  Enhancement  methods  for  the  Yagi-like  antennas  printed  in  a 
multi-layer  anisotropic  medium. 
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-  Extension  of  the  Evolutionary  Programming  technique  of  section  ID  for  gain  optimization  of  a 
Yagi-like  coaxial-fed  patch  array. 

-  Design,  fabrication  and  measurement  of  an  optimized  Yagi-like  coaxial-fed  patch  array. 

-  Investigation  of  other  feed  mechanisms  for  exciting  the  driver  element,  i.e.,  microstrip 
line,  proximity-coupled  line,  etc. 

In  addition,  a  number  of  research  tasks  are  planned  for  the  topic-2  project,  “New  design 
concepts  for  low-profile  dual-band  CP  microstrip  antenna  elements”.  Those  tasks  include: 

-Development  of  design  guidelines  for  the  proposed  dual  frequency  circularly-polarized 
electromagnetically  fed  microstrip  antenna  configurations. 

-  Characterization  of  the  radiation  parameters  of  each  proposed  configuration. 

-  Fabrication  and  experimental  verification  of  selected  design  for  operation  in  GPS  and/or 
MSAT  frequency  bands. 
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APPENDIX  A 


The  mixed-potential  Green's  functions  of  a  multi-layered  microstrip  structure,  shown  in 
Figure  A.l,  are  of  the  Sommerfeld  integral  type  which  can  be  expressed  as  the  inverse  Bessel 
transformation, 

OO 

G  Jz,p)  =  2n  fGg'Jz.X)J0(k(p\)hlk  ;  p  =  ^(x-x')2  +  (y-y')2  (A.l) 

o 


where  the  spectral-domain  Green’s  functions  G  H  (z,  A)  satisfy  the  voltage  distribution  in  the 
equivalent  multi-section  transmission  line  of  Figure  A.2, 


GEi  -  —2  (Vmi-Vei)  ,‘  GMi  -  Vmi 

A 

The  characteristic  impedance  in  each  section  is  given  by: 

Zei  =  (Z4  =  —  ,  for  Ve  ;  Zm,  =  (ZTE),  =  ^ 

(OEri  y  i 

y  i  =  -;A0V A  -  £  ri  M-  ri 


(A.2) 


for  Vm 


where  Yi  is  the  propagation  constant  in  the  i-th  section  of  the  line.  In  Figure  A.2,  Zgs  is  the  surface 
impedance  of  the  ground-plane;  for  an  structure  without  a  ground-plane,  Zgs  should  be  replaced  by 
an  open  circuit  representing  the  semi-infinite  space. 

To  derive  semi-analytical  solutions  for  G  we  first  extract  out  the  static,  1/  kop  singularity  of 
G;  the  remaining  integral,  gr ,  is  then  divided  into  four  regions  as  shown  in  Figure  A.3:  real-axis 
regions  from  0  to  Xs  (=  1-  8i)>  ^e  (=  nmax-  S2)  to  some  kmax  ( »  1)  and  Xmax  to  infinity,  and  the 
complex-contour  region  between  ks  and  ke.  Here  nmax  is  to  the  largest  refractive  index  of  the 
dielectric  layers,  and  8j  and  82  are  chosen  as  some  arbitrary  small,  but  not  very  close  to  zero, 
positive  values.  We  note  that  this  contour  avoids  any  possible  surface-wave  poles,  located  on  the 
real  axis  in  the  region  1  <  A  <  nmax .  Cubic-spline  interpolations  are  used  to  curve-fit  the  integrand, 
G,  in  each  region,  except  in  the  Amax  to  infinity  region  for  which  the  integration  range  is  divided 

into  small  pieces,  using  an  adaptive  marching  scheme,  and  G  is  approximated  on  each  piece  by  the 
first  few  terms  in  it’s  Taylor-series  expansion  in  A"^.  Care  should  be  exercised,  however,  in 
selecting  the  height  of  the  deformed  contour,  h,  above  the  surface-wave  poles  in  order  to  ensure 
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the  smoothness  of  the  integrand,  G,  in  (A.l)  before  applying  the  cubic-spline  interpolation. 
Numerical  experimentation  of  various  structures  show  that  accurate  results  for  the  deformed  part  of 
the  contour  can  be  obtained  by  selecting  the  height  h  in  Figure  A.3  as  Min(0.2,  2/kopmax)  <  h  < 

Max(0.8, 2/kQpmax)  in  the  cubic-spline  interpolation. 

Once  the  curve-fit  as  well  as  the  Taylor-series  coefficients  are  determined,  the  resulting 
integrals  of  the  type, 

Qm  =  fXm  J0(k0pX)dX  ,  m  =  -  4,  -  2,  0,  1  ,  2 , 3 .  (A.3) 

can  be  expressed  in  closed  forms  in  terms  of  Bessel  and  Struve  functions.  For  the  first  two  real- 
axis  regions,  the  final  expression  for  gr(p)  is  of  the  form  [14], 


,  ,  v/[(A.- An-/)  ,  . 

gr(p)= Xl — ; - 2o(‘^x») 

tTl[  koP 


( Bn  —  Bn- 1  )  _  .  ( Cn  — Cn  -1  )  _  /  (Dn  ~Dn  -1  )  -i 

+ - ^ — *2+;(VM+  — - — Q+2  (Vx«)  -  77  Q+i 


m3 


(Kp) 


(A.4) 


where  An,  Bn  and  Cn  are  functions  of  the  cubic-spline  coefficients,  and 

Q0(x)  =  1-xJ0(x)  +  ?y[J0(x))\i(x)-J1(x)Y\(J(x)}  ;  Q+I(x)  =  xJ,(x) 

Q+2(x)  =  Q0(x)  +  xJ0(x)  +  x2Jj(x)  ;  Q+3  (x)  =  2x2J0(x)-x(x2  -4)J,(x)  (A-5) 

wherein  Jj  and  Hj  are  the  Bessel  and  Struve  functions  [13],  respectively.  The  expression  for  the 
Xmax  to  infinity  region  is  similar  to  (A.4)  but  in  terms  of  Qo,  Q-2  and  Q-4  [14]. 

We  note  that  Qm  expressions  are  functions  of  distance  p  only  and  independent  of  the 
dielectric  layers'  thickness  and  material  properties.  Therefore,  once  the  constants  in  the  cubic- 
spline  interpolations  and  Taylor-series  expansions  are  computed  and  stored,  the  Sommerfeld 
integral  can  be  evaluated  for  all  values  of  p  at  once,  resulting  in  a  substantial  reduction  in  the  CPU 
time.  Additional  CPU  time  saving  can  also  be  achieved  for  the  first  three  regions  of  the  contour  if 
one  notes  that  the  Qm  functions  in  (4)  need  to  be  computed  only  once  for  a  given  set  of  kop 
values,  stored  and  then  used  for  all  M  Green’s  functions  of  an  N-layer  structure. 
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